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Abstract. This paper deals with mathematical issues relating to the computation 
of spectra of self adjoint operators on Hilbert spaces. We describe a general method 
for approximating the spectrum of an operator A using the eigenvalues of large finite 
dimensional truncations of A. The results of several papers are summarized which 
imply that the method is effective in most cases of interest. Special attention is paid 
to the Schrodinger operators of one-dimensional quantum systems. 
\ We believe that these results serve to make a broader point, namely that numerical 

problems involving infinite dimensional operators require a reformulation in terms of 
C*-algebras. Indeed, it is only when the given operator A is viewed as an element 
of an appropriate C* -algebra A that one can see the precise nature of the limit of 
the finite dimensional eigenvalue distributions: the limit is associated with a tracial 
state on A. For example, in the case where A is the discretized Schrodinger operator 
associated with a one-dimensional quantum system, A is a simple C*-algebra having 
I" ■ a unique tracial state. In these cases there is a precise asymptotic result. 
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1. Introduction. We discuss methods for computing the spectrum a (A), and 
especially the essential spectrum cr e (A), of bounded self-adjoint operators A on 
separable Hilbert spaces. 

We mean to take the term 'compute' seriously. Most operators that arise in 
practice are not presented in a representation in which they are diagonalized, and 
it is often very hard to locate even a single point in the spectrum of the operator 
[6], [7], [8], [12], [13], [15]. Some typical examples will be discussed in section 3. 
Thus, one often has to settle for numerical approximations to cr(A) or a e (A), and 
this raises the question of how to implement the methods of finite dimensional 
numerical linear algebra to compute the spectra of infinite dimensional operators. 
Unfortunately, there is a dearth of literature on this basic problem and, so far as 
we have been able to tell, there are no proven techniques. 

In this paper we establish an effective method for approaching such problems; we 
discuss issues associated with operator theory and operator algebras, but not issues 
belonging properly to numerical analysis. Thus, we address the question of whether 
or not certain finite dimensional approximations converge to the correct limit, but 
we do not address questions relating to how fast they converge nor how this method 
might be implemented algorithmically. Nevertheless, it may be appropriate to point 
out that these methods have been effectively implemented in a Macintosh program 
[2] which is available from the author. 

The following section contains a general method for computing the essential spec- 
trum of a self adjoint operator. The method depends on choosing an orthonormal 
basis and approximating the operator with n x n matrices obtained as sections of 
the infinite matrix of the operator with respect to this basis. Not every orthonor- 
mal basis is appropriate, and we develop criteria which show how the basis must be 
chosen. This involves an abstraction of the classical notion of band-limited matrix. 

Our method is effective for computing the essential spectrum. When there is a 
difference between the spectrum and essential spectrum, our results give no infor- 
mation about points in the difference a (A) \ a e (A). However, since a (A) \ a e (A) 
consists merely of isolated eigenvalues of finite multiplicity, this difference is usually 
insignificant and often empty (as it is for the main examples below). 

In section 3 we describe the source of our principal examples of self adjoint 
operators. We show how one should 'discretize' the canonical commutation relations 
so as to preserve the uncertainty principle, and relate the resulting Schrodinger 
operators to the tridiagonal operators considered in section 4. The one-parameter 
family of C*-algebras generated by the discretized canonical commutation relations 
turns out to be the one-parameter family of noncommutative spheres of Bratteli, 
Elliott, Evans and Kishimoto [9], [10], [11]. In our case the parameter is related to 
the numerical step size. 

In section 4 we clarify the role of C*-algebras in numerical problems of this 
kind. In particular, the C*-algebra s associated with a broad class of tridiagonal 
operators are simple C*-algebra s having a unique tracial state. The tracial state 
plays an essential role in describing the limit of the eigenvalue distributions of the 
approximating sequence of finite dimensional truncations of the basic operator. 

While there are some new results below, what follows is primarily an exposition 
of the results of several recent papers [1] , [2] , [3] , [5] . 
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increasing sequence 

T={H 1 QH 2 Q...} 
of finite dimensional subspaces H n of H with the property that 

lim dimH n = oo. 

n^oo 

The filtration T is called proper or improper according as the union UH n is dense 
in H or is not dense in H. In general, we will write P n and P+ for the projections 
onto the subspaces H n and UH n respectively. 

The simplest nitrations are associated with orthonormal sets in H. For example, 
if {e n : n = 1, 2, . . . } is an orthonormal set then 

(2.1) H n = [e 1 ,e2,...,e n ], n=l,2, ... 

defines a filtration of H which is proper iff {e n } is an orthonormal basis. If, on the 
other hand, we are given a bilateral orthonormal basis {e n : n = 0, ±1, ±2, . . . } for 
H, then 

H n = [e n , e_ n -|_i, . . . , 6 n ] 

defines a proper filtration in which the dimensions increase in jumps of 2. Moreover, 
there is a natural improper filtration associated with such a bilateral basis, namely 
F+ = {H+ CH+C...} where 

H n = [e ,ei,...,e n ], n = 1,2, .... 

The spaces ii^~ span a proper subspace -H+ of H. This example of an improper 
filtration is associated with 'unilateral' sections. That is to say, the matrix of any 
operator A e B(H) relative to the basis {e n : n G Z} is a doubly infinite matrix 
(a^j), whereas the matrix of the compression 

P+A \ H+ 

is a singly infinite submatrix of (a^ ). While these are the main examples of nitra- 
tions for our purposes here, nitrations in general are allowed to have very irregular 
jumps in dimension. 

Given any filtration T = {H n : n > 1} of H, we introduce the notion of degree 
of an operator (relative to J 7 ) as follows. 

Definition 2.2. The degree of A e B(H) is defined by 

deg(A) = sup rank(P n A — AP n ). 

n 

The degree of an operator can be any nonnegative integer or +oo, and it is clear 
that deg(A) = deg(A*). Operators of finite degree are an abstraction of band- 
limited matrices. For example, suppose we have a proper filtration arising from an 
orthonormal basis as in 2.1. If the matrix (a^) of an operator A e B(H) relative 
to this basis satisfies 
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then deg( A) < k. The degree function has a number of natural properties, the most 
notable being 

deg(AB) < deg(A) + deg(B). 

The set of all operators of finite degree is a unital *-subalgebra of B(H) which is 
dense in the strong operator topology [3]. 

There is a somewhat larger *-algebra of operators that one can associate to a 
filtration T . This is a Banach *-algebra relative to a new norm, and it plays a key 
role in the results to follow. This Banach algebra is defined as follows. Suppose 
that an operator A E B(H) can be decomposed into an infinite sum of finite degree 
operators 

(2.3) A = A 1 + A 2 + ... 
where the sequence A n satisfies the condition 

oo 

(2.4) ^^(l + deg^) 1 / 2 )!!^!! <oo. 

n=l 

Notice that the sum indicated in 2.3 is absolutely convergent in the operator norm 
since 2.4 implies 

oo 

||-A n || < S < OO. 

n=l 

We define \A\jr to be the infimum of all numbers s associated with decompositions 
of A of this kind. If A cannot be so decomposed then \A\jr is defined as +oo. We 
define 

V T = {A E 13(H) : \A\jr < oo}. 

(7>r, | • | jit) is a unital Banach algebra which contains all finite degree operators, 
and the operator adjoint defines an isometric involution in V^. The two norms 
currently in view are related by 

\\M < \A\ T . 

While it is not so easy to calculate the norm in "Djf, it is usually very easy to 
obtain effective estimates. For example, let T be a proper filtration associated with 
a unilateral orthonormal basis as in 2.1, let A be an operator on H and let (a^) be 
the matrix of A relative to this basis. Letting 

d k = sup |a M+fe | 

i 

denote the supnorm of the kth diagonal, k = 0, ±1, ±2, . . . , then we have 

+ oo 

\A\ T < (1 + (2|^|) 1/2 K. 

k= — oo 

See [3]. Thus, A will belong to D^r whenever the diagonals of (a^) die out fast 
enough so that 

+oo 

\k\ 1/2 d k < oo. 
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We now indicate the role of the Banach algebra Dp in computing spectra. Let 
A E B(H) be a self adjoint operator and consider the sequence of matrices A n 
obtained by compressing A along the filtration T = {Hi C H 2 C . . . }: 

A n = P n A \ Hn . 

We are interested in certain asymptotic quantities that can be computed (at least 
in principle) from the sequence of finite dimensional spectra a(Ai), (7(^2), .... The 
simplest one is the set A, which consists of all real numbers with the property that 
there is a sequence of eigenvalues A n G a(A n ), n = 1, 2, . . . such that 

lim X n = A. 

n^oo 

A is in some sense the set of limit points of the sequence of sets a(Ai), a(A 2 ), .... 
Notice however that A is smaller than the topological limit superior of sets because 
limits along subsequences of a(Ai), (7(^2), ... do not qualify for membership in A. 
It is easy to see that A is closed, and in general A contains the spectrum of A [3]. 
In particular, A is never empty. The points of A are classified as follows. 

For every subset SCI and every n = 1, 2, . . . we write N n (S) for the number of 
eigenvalues of A n which belong to 5\..and of course one counts multiple eigenvalues 
according to their multiplicity. 

Definition 2.5. A real number A is called a transient point if there is an open set 
U containing A such that 

sup N n (U) < M < 00. 

n>l 

A is called an essential point if for every open set U containing A we have 

lim N n (U) = 00. 

n^oo 

Remarks. With every transient point A of A there is an associated pair of positive 
integers p < q with the property that for every sufficiently small neighborhood U 
of A one has the behavior 

P<N k (U) < q 

for large enough k > n = njj, and moreover both extreme values p and q are taken 
on infinitely many times by the sequence Ni(U), AT 2 (t/), .... 

The set of all essential points is a subset of A which we denote by A e . Again, 
it can be seen that A e is a closed set which contains the essential spectrum of A 
[3]. At this level of generality, there does not appear to be much more that one 
can say. For instance, there are examples which show that the both inclusions 
cr(A) C A and a e (A) C A e can be proper (see the appendix of [3]). Other examples 
show that A can contain points which are neither essential nor transient; for every 
small neighborhood U of such a point one can find subsequences n\ < n 2 < ... 
for which the sequence of positive integers A^ ni (U), N n2 (U), . . . stays bounded, and 
other subsequences mi < ni2 ■ ■ ■ such that N mk (U) — > 00 as k — > 00. Fortunately, 
the following result implies that in the reasonable cases we will not find this kind 
of instability. 

We can now state our main general result, which shows how one must choose a 
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Theorem 2.6. Let T be a proper filtration and let A be a self adjoint operator 
which belongs to the Banach *-algebra V^. Then A e coincides with the essential 
spectrum of A. Moreover, every point of A is either transient or essential. 

Remarks. In the following section we will encounter tridiagonal operators A which 
are defined in terms of a bilateral orthonormal basis {e n : n = 0, ±1, ±2, . . . } by 

(2.7) Ae n = e n _i + d n e n + e n+1 , 

where d n represents a bounded sequence of reals. Theorem 2.6 shows how the 
essential spectrum of A can be computed in terms of the eigenvalue distributions 
of the sequence of (2n + 1) x (2n + 1) matrices obtained by compressing A along 
the sequence of subspaces H n = [e_ n , e_ n+ i, . . . , e n _i, e n ], n = 1,2,.... While 
Theorem 2.6 says nothing about rates of convergence, our experience with the 
operators of section 4 has shown that convergence is rapid; in fact, it is fast enough 
to allow the construction of excellent pictures of cr e (A) on desktop computers in a 
minute or two [4]. Nevertheless, there remains an important problem of obtaining 
an appropriate definition of the "rate" of convergence of such approximations, and 
the estimation of this rate. The paper [3] does not address these issues of error 
definition and estimation. 

In carrying out computations, it is usually more convenient to deal not with 
the above "bilateral" sections A n but with smaller "unilateral" sections. More 
precisely, these are defined in terms of an improper filtration {H+ C H% C . . . }, 
= [ei, e2, • • • , e n ], and the corresponding compressions of the operator A 

A n = \h+' 

.P+ denoting the projection onto H+. Since the filtration is improper, Theo- 

rem 2.6 does not allow one to draw conclusions about the essential spectrum of A, 
but rather the essential spectrum of the operator A + obtained by compressing A 
to the subspace H + generated by \J n H+. Thus there remains a significant problem 
of relating the essential spectrum of A to that of A + . These issues will be taken up 
in section 4 below. 

3. Discretized Schrodinger operators. In this section we will indicate how 
a one-dimensional quantum system should be discretized in order to carry out 
numerical computations. We will find that the "discretized" canonical commutation 
relations generate a C*-algebra which is isomorphic to one of the noncommutative 
spheres of Bratteli, Elliott, Evans and Kishimoto (with parameter related to the 
numerical step size), and we will find that the resulting discretized Hamiltonian is 
a bounded self-adjoint operator which is amenable to the methods of the preceding 
section. 

Most one-dimensional quantum systems are modelled on the Hilbert space L 2 (IR). 
The canonical operators are the unbounded self-adjoint operators P, Q defined on 
appropriate domains in L 2 (M) by 

1 d 

i dx ' 



(3.1) 



P 
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They obey the canonical commutation relations on an appropriate common domain 

PQ-QP=-1. 

i 

The time development of the system is described by a one-parameter unitary group 
of the form 

W t = e UH , t E K, 
where H is the Hamiltonian of the system 

(3.2) H=±P 2 + <I>(Q), 

4> being a real-valued continuous function of a real variable which represents the 
potential of the classical system being quantized. 

In order to carry out numerical computations one first has to replace the differ- 
ential operator H with an appropriately discretized version of itself. Moreover, one 
has to decide how this should be done so as to conform with the basic principles 
of numerical analysis while at the same time preserving the essential features of 
quantum mechanics (i.e., the uncertainty principle). In [1] we presented arguments 
which we believe justify the following procedure. 

One first settles on a numerical step size a. This can be regarded as a small pos- 
itive rational number, whose size represents the smallest time increment to be used 
in the difference equations that replace differential equations. One then discretizes 
both operators P and Q, and finally uses the formula 3.2 to define the corresponding 
discretized version of the Hamiltonian. 

In more detail, we replace the differential operator P with the bounded self- 
adjoint difference operator P a defined by 

Zia 

Noting that the one-parameter group of translations Vtf(x) = f(x + t) is generated 
by P in the sense that Vj = e ltp , t G K, we have 

P a = -^-(e iaP - e~ iaP ) = - sin(aP). 

Now we must discretize Q but we must be careful to do it in a way that preserves 
the uncertainty principle insofar as that is possible. In section 3 of [1], we argued 
that this requirement imposes a very strong restriction on the possible choices for 
the discretized Q, and that in fact the only "correct" choice is given by 

Qa = ^-(e 1 ^ - e~ i(J Q) = - sin(aQ). 

More explicitly, we have 

Q_f( x ) = -p.m(ax)f(x]. x(=U.f(= L 2 (M). 
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The resulting discretized Hamiltonian is then defined as follows: 



(3-3) H a = -P 2 + V(Q a ). 



Once we have the operator H a we are in position to carry out numerical com- 
putations with the quantum system. For example, if the system is described at 
time t by a wave function / e L 2 (R) then the state of the system at time t + At is 
approximated by the wave function 



We caution the reader that, while the preceding formula is convenient for illus- 
trating one way to make use of the discretized Hamiltonian, it must be modified 
appropriately in order to correctly model the dymanical group in practice because 
the operator / i— > g is not unitary. Readers interested in carrying out numerical 
computations can find a discussion of closely related issues in [14, pp 662-663]. 

Let T> a be the C*-algebra generated by the set of operators {Per, Qa}- T) a is the 
norm-closed linear span of all finite products of terms involving either P a or Q a . It 
is not obvious that V a contains the identity operator but that is the case. V a is the 
discretized counterpart of the algebra of observables, and notice that it contains 
the operator H a . Thus it is important to understand the structure of T> a . 

In fact, T> a is a simple unital C*-algebra which is isomorphic to one of the 

non-commutative spheres of Bratteli et al [9, 10, 11]. Moreover, while the oper- 
ators Per, Q a no longer satisfy the canonical commutation relations, they do obey a 
more subtle discretized form of the CCRs, and more generally it is the universal C*- 
algebra associated with these "discretized" CCRs that is naturally associated with 
the non-commutative spheres. The reader is referred to [2] for a detailed discussion 
of these and related issues. 

In particular, if one is interested in computing the spectrum of H a then one is 
free to choose any convenient representation of V a and compute the spectrum of 
H a in that representation. Since V a is simple, the spectrum does not depend on 
the representation chosen. Actually, the most convenient realization of H a is one 
in which it is a tridiagonal operator. In this case, it is more appropriate to work 
with a subalgebra of V a which contains H a . The precise statement follows. 

Proposition 3.4. Let A be the C* -algebra generated by P 2 and Q a , and let K be 

a Hilbert space spanned by a bilateral orthonormal set {e n : n e Z}. Then there is 
a faithful representation it : A — > B(K) such that 7r(H a ) has the form 



g=(l + iAtH a )f = f + iAtH a f. 



(3.5) 



7r(H a ) = aT + bl 



where a = l/8o~ 




and T is the tridiagonal operator 



(3.6) 



Te n = 



e n _i + 8a 2 (f)(- sin(2ncx))e n + e n+1 , 



a 



n = 0,±l,±2, .... 

proof. Consider the unitary operators U, V defined by 



Uf(x) = e™f{x) 



C*- ALGEBRAS IN NUMERICAL LINEAR ALGEBRA 



9 



Then 



and notice that 



Q a = ^-(U-U-^ 



Pl = - _L (yi + V- 2 )+ ' 



4a 2 v ' 2a 2 ' 

Let B be the C*-algebra generated by U and V 2 . Clearly A is contained in B, and 
because of the commutation relation 



V 2 U = e 2ia2 UV 2 



and the fact that a is a positive rational number, it follows that B is an irrational 
rotation C*-algebra. We will define a representation tt\ of B on K; the required 
representation of A is obtained by restriction. 

In order to specify a representation of B on K, it is sufficient to specify a pair 
of unitary operators S, D on K satisfying 



(3.7) SD = e 2ia2 DS. 

7ri is then uniquely defined by specifying that 7ii(V 2 ) = S, ivi(U) = D. Let 



2 



n = 0, ±1, ±2, .... Clearly £ and D are unitary operators and the reader can verify 
3.7 directly. Hence there is a representation 7Ti of £> with the stated properties. 
Noting that 

1 s -JL s-i + JL, 



4(T 2 4(T 2 2(T 2 

and that 7ri(Qo-) is the diagonal operator 

7ri(Q CT )e n = -^(e 2 ^ 2 - e - W )e n = - sin(2na 2 )e n , 
2zcr a 

we may conclude that 

^i((/}(Qa))e n = (f>(- sin(2na 2 ))e n 

(7 

for every n = 0, ±1, ±2, .... Combining these two formulas, we find that the image 
of H a is given by 
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which has the form spelled out in 3.5 and 3.6 □ 
Remark 3.8. Consider the diagonal sequence 

d n = 8a 2 (p(- sin(2na 2 )), n G Z 
a 

appearing in the formula 3.6. We want to point out that if the function is 
continuous and not constant on the interval [—-,+-1 and if a 2 is not a rational 
multiple of n, then the sequence (d n ) is almost periodic but not periodic. Indeed, 
the sequence n t— > sin(2na 2 ) is almost periodic because it is a linear combination 
of complex exponentials of the form e ian where a is a real number. Since the set 
of all almost periodic sequences form a commutative C*-algebra it is closed under 
the continuous functional calculus, and therefore the sequence (d n ) must be almost 
periodic. Note too that (d n ) cannot be periodic. For if there did exist integers 
p > 1 and n such that 

d-n+kp = d n , k = 0, ±1, ±2, . . . 

then since a 2 is not a rational multiple of it the numbers 

sin(2na 2 + Ikpa 2 ), keZ 

would fill out a dense set in the interval [—1,-1-1], and hence 4> would have to be 
constant on the interval [— -, 

The material presented in this section leads toward a significant conclusion. The 

problem of computing the spectrum of the discretized Hamiltonian of a one dimen- 
sional quantum system can be reduced to the problem of computing the spectrum of 
a self-adjoint tridiagonal operator of the form 

(3.9) Te n = e n _i 

~\~ d n e n + e n+ i, nGZ 

where {d n : n G Z} is a bounded almost periodic sequence of reals which is not 
periodic. For such operators one can work with either unilateral sections or bilateral 
sections, as we will see in the following section. For example, if for n > 1 we let T n 
be the compression of T onto the linear span of {ei, e2, . . . , e n } then even though we 
are in effect working with an improper filtration, one may apply a suitable variation 
of Theorem 2.6 to conclude that the spectrum of T is the set of all essential points 
associated with the sequence of self-adjoint matrices Ti, T 2 , . . . . 

In the following section, we will obtain more precise information about the dis- 
tribution of the eigenvalues of the sequence T n , n = 1,2, 

4. Limits, simple C*-algebras, and traces. Let (d n ) ne z be a bounded almost 
periodic sequence of reals which is not periodic and let T be the tridiagonal operator 
of 3.9 

(4.1) Te n = e n _! -\- d n e n + e n+1 , n E Z, 

{e n : n G Z} being an orthonormal basis for a Hilbert space H. We will show that 
the eigenvalue distributions of the n x n sections of T actually converge (in the 
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becomes important to understand the nature of this limiting measure \it- We show 
that \it is associated with a tracial state on a certain C*-algebra At associated with 
T. At certainly contains T but it is much larger than the C*-algebra generated by 
T; indeed, At is a simple C*-algebra having a unique tracial state. 
Let S be the bilateral shift defined on H by 



and let D be the diagonal operator associated with the sequence (d n ), 

De n = d n e n , n G Z. 

At is defined as the C*-algebra generated by both operators D and S. At is a 
separable unital C*-algebra which contains T, and of course At depends on the 
particular choice of diagonal sequence (d n ). When T is a discretized Hamiltonian 
as in the previous section, At will depend on both the numerical step size a and 
the potential <j>. Nevertheless, in all cases we have 

Theorem 4.2. At is a simple unital C* -algebra- subalgebra of 13(H) which has a 
unique tracial state. 

Remark. A tracial state of At is a linear functional r : At — > C satisfying 



for every 1,7 6 At, and which is normalized so that r(l) = 1. Theorem 4.2 is 
proved in proposition 3.2 of [5]. 

More generally, suppose we are given an arbitrary concrete C*-algebra A C 
13(H). We need to single out the nitrations that are "compatible" with A. Let 
T = {Hi C H 2 C . . . } be a filtration of H which may be improper, and put 



As in section 2 we may speak of the degree of an operator X e 13(H) relative to 
the filtration T. Following [3], we say that T is an ^4-filtration if the set of finite 
degree operators which belong to A is norm-dense in A. 

In connection with the operators T of 4.1, we will consider the filtration JF T = 
{H n } where H n = [e 1 , 

e 2? • • • j e n]- This is an improper filtration for which 



More significantly, we have 

Proposition 4.3. JF T is an At -filtration. 

4.3 is a simple consequence of the fact that the finite degree operators in At 
form a of At, and that the operators D and S have respective degrees and 1 (see 
the proof of Theorem 3.4 of [5]). For every n > 1 let T n be the compression of T 
to H n . Relative to the obvious basis for H n: the matrix of T n is 




n e Z 



t(X*X) > 0, and 
t(XY) = t(YX) 



H+ = U n if, 



H+ = [ei,e 2 ,...]. 



/d 1 
1 




1 



1 





1 
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There are two issues that need to be understood. The first has to do with the 
relation between operators in At and their compressions to H+; the second requires 
relating the trace on At to the limit of the eigenvalue distributions of the sequence 
of matrices Ti, T2, . . . . 

In order to discuss the first of these, we consider a more general setting in 
which we are given a unital C*-algebra A C 13(H) and an improper filtration 
{Hi C H 2 C . . . }. We will consider the space A + of all operators on H + having 
the form 

P+A \ H+ +K, 

where A e A and K is a compact operator on H+, P+ denoting the projection on 
H + . We will write /C+ for the algebra of all compact operators on H + . 

Theorem 4.4. Assume that A has a unique tracial state r, let {Hi C H2 C . . . } 
6e an A-filtration and assume that H + has the following property 

(4.5) A \h + = compact =>- A = compact 

for every operator A e A. Then A + is a C* -algebra containing 1C+ which has a 
unique tracial state r+. t+ is related to r by way of 

r+(P+A \ H+ +K) = t(A), AeA,KelC+. 

Moreover, the natural map of A to the quotient A+ / JC + given by 

A ^ P+A \ H+ +/C+ 

is an isomorphism of C* -algebras : 

A = A+/K+. 

Remarks. The argument required here can be found in the proof of Theorem 2.3 
of [5]. These results depend on the following relationship that exists between the 
operators in A and the projection P + associated with an ^4-filtration: 

A <E A =>- P+A - AP+ G K 

see [5, proposition 2.1]. 

Theorem 4.4 implies that we have a short exact sequence of C*-algebras 

-> /C+ -> A+ -> ^ -> 0. 

This shows that the structure of is somewhat analogous to the structure of the 
Toeplitz C* -algebra T, 

-> /C -> T -> C(T) -> 0, 

except that in our applications the quotient C*-algebras At are highly noncommu- 
tative. 

Finally, we remark that the hypothesis 4.5 is satisfied for our examples At 
because At is a simple C*-algebra with unit and H+ is infinite dimensional (see 
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Now we can define the measure \xt alluded to at the beginning of this section. 
Let r be the unique tracial state of At- Using the functional calculus for bounded 
self-adjoint operators, we can define a positive linear functional on Cq(R) by 



feC 



r(f(T)). 



By the Riesz- Markov theorem, there is a unique positive measure \it on R such 
that 



/+oo 
f(x)diir(x) = T(f(T)), feC 
-oo 



Ht is obviously a probability measure whose support is contained in the spectrum 
of T. ht is called the spectral distribution of T. Because At is simple r must be a 
faithful trace, and hence the closed support of ht is exactly the spectrum of T. 

The preceding discussion, together with the general results of [3, section 4] can 
be applied to obtain the following result, which is Theorem 3.4 of [5]. 

Theorem 4.6. For every positive integer n let A™ < X n < • • • < A™ be the eigen- 
value list of the symmetric n x n matrix 



(4.7) 



Then for every f e Cq(K) we have 





1 





1 


d 2 


1 


; 


1 


d 3 











Vo 








.) we 


have 






















• • • d n -i 


1 


1 


d n ' 



lim -(/(A?) + /(A?) + ...,+/(AS)) 

n— >oo 77, 



+oo 



f(x)d/i T (x). 



Remarks. Theorems 4.6 and 2.6 together provide rather precise information about 
the rate at which the eigenvalues of the matrices 4.7 accumulate at points in and 
out of the spectrum of T. For example, let A G cr(T) and let / be an open interval 
containing A. Then /j,t(I) > 0, and if a and (3 are chosen close to /j,t(I) in such a 
way that a < Ht(I) < then the number N n (I) of eigenvalues of T n which belong 
to / will satisfy the inequalities 

an < N n (I) < (3n 

for all sufficiently large n. 

For the operators T of 4.1, it is not hard to show that the essential spectrum 
of T is identical with cr(T). Thus we can apply Theorem 2.6 above to conclude 
that if A does not belong to cr(T) then for every sufficiently small open interval 
containing A, the sequence of numbers Ni(I), N 2 (I), . . . actually remains bounded. 
This behavior is quite visible in the pictures generated by the program [4]. 

Some additional remarks relating to computational issues can be found in the 
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